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Abstract 

The theoretical description of the correlations between observables in two separated 
rapidity intervals for AA-interactions at high energies is presented. In the case with 
the real nucleon distribution density of colliding nuclei the MC calculations of the 
long-range correlation functions at different values of impact parameter are done. 
For n-n and pt~n correlations it is shown that the impact parameter fluctuations 
■ at a level of a few fermi, unavoidable in the experiment, significantly change the 
magnitude of correlation coefficients. The rise of pt~n and especially pt~Pt correlation 
coefficients is found when one passes from SPS to RHIC and LHC energies. 
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1 String fusion model (SFM) 

in 

The colour string model |IJ Ej originating from Gribov-Regge approach is being widely 
applied for the description of the soft part of the multiparticle production in hadronic and 
nuclear interactions at high energies. In this model at first stage of hadronic interaction 
Ph! the formation of the extended objects - the quark-gluon strings - takes place. At second 
stage the hadronization of these strings produces the observed hadrons. In the original 
version the strings evolve independently and the observed spectra are just the sum of 
individual string spectra. However in the case of nuclear collision, with growing energy 
and atomic number of colliding nuclei, the number of strings grows and one has to take 
into account the interaction between them. 

One of possible approaches to the problem is the colour string fusion model The 
model is based on a simple observation that due to final transverse dimensions of strings 
they inevitably have to start to overlap with the rise of their density in transverse plane. 
At that the interaction of string colour fields takes place, which changes the process of 
their fragmentation into hadrons as compared with the fragmentation of independent 
strings. So we have one more interesting nonlinear phenomenon in nuclear interactions 
at high energies - the field of physics the investigations in which were initiated by pioneer 
works of academician A.M. Baldin 

It was shown |H] that the string fusion phenomenon considerably damps the 
charged particle multiplicity and simultaneously increase their mean p t value as compared 
with the case of independent strings. In accordance with a general Schwinger idea [7] and 
the following papers [HI E] (colour ropes model) two possible versions of string fusion 
mechanism were suggested. 

The first version [S] of the model assumes that the colour fields are summing up only 
locally in the area of overlaps of strings in the transverse plane. So we will refer to this 
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local fusion or overlaps. In this case one has 

(n) k = ^Vk (p 2 t ) k = P 2 Vk (1) 

e>"o 

Here (n) k is the average multiplicity of charged particles originated from the area S k , where 
k strings are overlapping, and (p 2 ) k is the same for their squared transverse momentum. 
The Ho and p 2 are the average multiplicity and squared transverse momentum of charged 
particles produced from a decay of one single string, and <To is its transverse area. 

In the second version PU| of the model one assumes that the colour fields are summing 
up globally - over total area of each cluster in the transverse plane - into one average 
colour field. This case corresponds to the summing of the source colour charges. We will 
refer to this case as a global fusion or clusters. In this case we have 

(n) d = HQ—\Jkd (p 2 t ) cl = p 2 \J~k~ d k d = Nd a ° (2) 

Here (n) d is the average multiplicity of charged particles originated from the cluster of 
the area S d and (p 2 ) d is the same for their squared transverse momentum. The N d r is 
the number of strings forming the cluster. 

Note that in two limit cases both versions give the same results. 

For N non-overlapping strings we have in the local version: k = 1, Si = Nao, (n) = 
(n) 1 = Nfi and (p 2 ) = (p 2 )i = p 2 ■ In the global version in this case we have N clusters 
each formed by only one string, so k d = 1, (n) = N(n) d = Nfi and (p 2 ) = {p 2 ) d = p 2 ■ 

For iV totally overlapped strings we have in the local version: k = N, Sn = ctq, 
(n) = (n) N = y/N /i and (pf) = (p 2 ) N = VN p 2 . In the global version in this case we 
have one cluster of the area S d = cr formed by N string, so k d = N, (n) = (n) d = VW /j 
tmd(p 2 ) = (p 2 t ) d = VNp 2 . 

So in both versions of the model when we pass from iV non-overlapping strings to iV 
totally overlapped strings the average multiplicity decreases from (n) = Nfi to (n) = 
\/~N and the mean p 2 t increases from (p 2 ) = p 2 to (p 2 ) = y/N p 2 . 



2 Cellular analog of SFM 

To simplify calculations in the case of real nucleus-nucleus collisions a simple cellular 
model originating from the string fusion model was proposed JT]. In the framework of 
the cellular analog along with the calculation simplifications the asymptotics of correlation 
coefficients at large and small string densities can be found analytically in the idealized 
case with the homogeneous string distribution, which enables to use these asymptotics 
later for the control of the Monte-Carlo (MC) algorithms. 

Two versions of the cellular model as the original SFM can be formulated - with local 
and global string fusion. In this model we divide all transverse (impact parameter) plane 
into sells of order of the transverse string size o~$. 

In the version with local fusion the assumption of the model is that if the number of 
strings belonging to the zj-th cell is kij, then they form higher colour string, which emits 
in average fj-oykij particles with mean p\ equal to p 2 ^Jk~ij (compare with (JT|)). Note that 
zero "occupation numbers" kij = are also admitted. 
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In the version with global fusion at first we define the neighbour cells as the cells with 
a common link. Then we define the cluster as the set of neighbour cells with non zero 
occupation numbers k^ 7^ 0. After that we can apply the same formulae of the global 
fusion (J2J as in the original SFM, where N^ r is the number of strings in the cluster and 
Set/a® is the number of cells in the cluster. 

From event to event the number of strings k^ in the ij-th cell will fluctuate around 
some average value - k^. Clear that in the case of real nuclear collisions these average 
values k^ will be different for different cells. They will depend on the position (s^) of the 
ij-th cell in the impact parameter plane (s is two dimensional vector in the transverse 
plane). In the case of nucleus-nucleus A5-collision at some fixed value of impact param- 
eter b one can find this average local density of primary strings k^ in the point using 
nuclear profile functions T^s^ + b/2) and 7b (s^- — b/2). 

In MC approach knowing the k^ one can generate some configuration C = {kij}. 
To get the physical answer for one given event (configuration C) we have to sum the 
contributions from different cells in accordance with local or global algorithm (see above), 
which corresponds to the integration over s in transverse plane. Then we have to sum 
over events (over different configurations C). Note that as the event-by-event fluctuations 
of the impact parameter at a level of a few fermi are inevitable in the experiment one has 
to include the impact parameter b into definition of configuration C = {b, k^}. 

3 Long-range correlations 

The idea EH Ej to use the study of long-range correlations in nuclear collisions for 
observation of the colour string fusion phenomenon based on the consideration that the 
quark-gluon string is an extended object which fragmentation gives the contribution to 
wide rapidity range. This can be an origin of the long-range correlations in rapidity 
space between observables in two different and separated rapidity intervals. Usually in 
an experiment they choose these two separated rapidity intervals in different hemispheres 
of the emission of secondary particles one in the forward and another in the backward in 
the center mass system. So sometimes these long-range rapidity correlations are referred 
as the forward-backward correlations (FBC). 

In principle one can study three types of such long-range correlations: 

n-n - the correlation between multiplicities of charged particles in these rapidity intervals, 

Pt-Pt - the correlation between transverse momenta in these intervals and 

p t -n - the correlation between the transverse momentum in one rapidity interval and the 
multiplicity of charged particles in another interval. 

Usually to describe these correlations numerically one studies the average value (B) F of 
one dynamical variable B in the backward rapidity window Ay B , as a function of another 
dynamical variable F in the forward rapidity window Ay F . Here (...) F denotes averaging 
over events having a fixed value of the variable F in the forward rapidity window. The 
(...) denotes averaging over all events. So we find the correlation function (B) F = f(F). 

It's naturally then to define the correlation coefficient as the response of (B) F on the 
variations of the variable F in the vicinity of its average value (F). At that useful also to 
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go to the relative variables, i.e. to measure a deviation of F from its average value (F) in 
units of (F), and the same for B. So it's reasonable to define a correlation coefficient b B _ F 
for correlation between observables B and F in backward and forward rapidity windows 
in the following way: 

L _ (F) d(B) F 



>B-F 



(B) dF 



(3) 



F=(F) 



As the dynamical variables we use the multiplicity of charged particles (n), produced 
in the given event in the given rapidity window, and the event(!) mean value of their 
transverse momentum (p t ), i.e. the sum of the transverse momentum magnitudes of all 
charged particles, produced in the given event in the given rapidity window (Ay), divided 
by the number of these particles (n): 



Pt = ~^2\Pti\ 



n 



where yi G Ay; 



i=i 



So we can define three correlation coefficients: 
(n F ) d(n B ) n 



(n B ) dn 1 



n F = { n F ) 



Pt —Pt 



= M_ d(p tB ) ni 
n ~ n ~ (PtB> dn F 



1, ...,n. 



(4) 



( PtF ) d( PtB ) Pti 



(PtB> d PtF 



(5) 

n F =i n F ) 

Here n B , n F are the multiplicities andp tB , p tF are the event (j3J) mean transverse momen- 
tum of the charged particles, produced in the given event correspondingly in the backward 
(Ay B ) and forward (Ay F ) rapidity windows. 



4 Results of the calculations 

In Figs. 1-3 the results of the MC calculations of these correlation coefficients are presented 
for nucleus-nucleus collisions at different values of the centrality. In all figures (o) and 
(•) denote the results of calculations in the framework of the original SFM (with the 
taking into account the real geometry of merging strings) for its local (overlaps) and global 
(clusters) versions correspondingly. The (□) and (■) denote the results of calculations in 
the framework of the cellular analog of SFM for its local and global (clusters) versions. 
All presented results are for the forward rapidity window of 2 unit length (Ay F = 2). The 
lines are only to guide the eye. 

In Fig.l we present the b n _ n correlation coefficient for AuAu collisions at RHIC energy 
and in Fig. 2 we present the b pt _ n correlation coefficient for PbPb collisions at LHC energy. 
In both figures the calculations are fulfilled three times: 

1) at fixed values of impact parameter (db = 0), 

2) with impact parameter fluctuations within 1 fm window (db = 1), 

3) with impact parameter fluctuations within the whole class of centrality (db = class) 
(for LHC by convention this value is taken to be equal 3 fm, db = 3) We see that the 
impact parameter fluctuations at a level of a few fermi significantly change the magnitude 
of correlation coefficients. Note also that all results obtained in the framework of the 
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Figure 1: The b n _ n correlation coefficient for AuAu collisions at ^/s = 130 GeV as a 
function of the impact parameter b for tree choices of impact parameter window db (see 
text). 

original SFM and its cellular analog for their local and global versions practically coincide, 
except for the p t -n correlation at fixed value of impact parameter at LHC energy, where 
the correlation coefficient b n is very small. 

In Fig. 3 the energy dependence of the b pt _ pt correlation coefficient is presented. The 
calculations are made for the 1 fm impact parameter window (db =1). In this case we 
see the considerable rise of p t -pt correlation coefficient from SPS to LHC energies. 

5 Conclusion 

In the case with the real nucleon distribution density of colliding nuclei the MC calcula- 
tions of the long-range correlation functions at different values of impact parameter are 
done. For n-n and Pt~n correlations it is shown that the impact parameter fluctuations at 
a level of a few fermi, unavoidable in the experiment, significantly change the magnitude 
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Figure 2: The same as in Fig.l but for the b pt _ n 
at = 5500 GeV. 



correlation coefficient for PbPb collisions 



of correlation coefficients for all centrality classes as compared to ones calculated earlier 
at the fixed values of impact parameter [T5j . 

It is shown also, that for the pt~Pt correlation the event-by-event correlation between 
event mean values of transverse momenta of the particles emitted in two different rapidity 
intervals does not decrease to zero with the increase of the number of strings in contrast 
with the correlation between the transverse momenta of single particles produced in these 
two rapidity windows which was studied earlier [T3] . 

The rise of pt~n and especially pt~Pt correlation coefficients is found when one passes 
from SPS to RHIC and LHC energies. 
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Figure 3: The b pt _ pt correlation coefficient as a function of the impact parameter b (at 
db = 1 fm) for tree choices of the initial energy: ^/s = 17; 130 and 5500 GeV. 
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